X(3872) is not a true molecule 
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The wave function of the charmonium-like meson X(3872) is expected to have a very significant 
D°D*° component, as the latter threshold lies only 0.16 MeV above the meson's mass. Some conclude 
from this mass coincidence that X(3872) is a meson-meson molecule, whatever the dynamics giving 
rise to the state. This would imply that the D°D*° component of the X(3872) wave function is the 
only relevant one. In the present paper we study this issue by employing a soluble model for a 3 Pi cc 
state coupled to an S-wave D°D decay channel, communicating via the s Po mechanism. The model 
is a simplified, coordinate-space version of the resonance-spectrum expansion previously employed 
in a detailed investigation of X(3872)'s nature. The resulting two-component wave function is 
calculated for different values of the binding energy (BE) and the transition radius a. Thus, a 
significant cc component is found in all situations. However, the long tail of the D°D*° wave 
function, in the case of small BEs, strongly limits the cc probability, which roughly lies in the range 
7-11%, for the true BE of 0.16 MeV and a between 2 and 3 GeV -1 . Furthermore, a reasonable 
value of 7.8 fm is obtained for the X(3872) r.m.s. radius at the physical mass, as well as an 5*-wave 
D° D scattering length of 11.6 fm. Finally, the S-matrix pole trajectories as a function of coupling 
constant show that X(3872) can be generated either as a dynamical pole or as one connected to the 
bare cc confinement spectrum, depending on details of the model. From these results we conclude 
that X(3872) is not a genuine meson-meson molecule, nor actually any other mesonic system with 
non-exotic quantum numbers, due to the inevitable mixing with the corresponding quark-antiquark 
states. 

PACS numbers: 14.40.Pq, 14.40.Lb, 12.39.Pn, 11.55.Ds 



I. INTRODUCTION 

The Belle Collaboration discovered [1] the charmon- 
ium-like state X(3872) almost a decade ago, observing it 
the decay B ± — > K ± tt + tt^ J/ip, with a significance in ex- 
cess of lOcr. The new meson was then confirmed by CDF 
0, DO 0, BABAR g], and recently also by LHCb [5]. 
In the 2012 PDG tables [BJ, A(3872) is listed with a mass 
of 3871.68 ± 0.17 MeV and a width T < 1.2 MeV. As for 
the quantum numbers, the possibilities are J PC = 1 ++ 
and J PC =2 h according to the PDG [BJ. The positive 
C-parity assignment resulted out from an analysis of an- 
gular distributions by CDF [7 . Recently, observation of 
the decay to 7J 'ftp by Belle [8[ unambiguosly confirmed 
C = +. For further experiments on A(3872) production, 
see Ref. [BJ. 

On the theory side, the discussion about the nature 
of A (3872) continues most vivid. Although the PDG |BJ 
lists 1 ++ and 2 h as the meson's possible quantum num- 
bers, and a recent BABAR analysis [5] of the 3ir invari- 



ant mass distribution in the decay X(3872) — > ujj/tp even 
seems to slightly favor the 2 h assignment, most model 
builders describe the state as an axial vector. For in- 
stance, model calculations of semi-inclusive B — s- f] C 2+X 
processes [TB] as well as electromagnetic rj c2 decays [TT] 
have been shown to disfavor the 2 h scenario. The same 
conclusion was reached in a tetraquark description of 
A (3872) [12], while pion exchange in a molecular picture 
would be repulsive in this case [13] and so inhibitive of 
a bound state. Finally, unqucnching a 1 x Di cc state by 
including meson-meson loops can only further lower the 
bare mass, which lies in the range 3.79-3.84 GeV for all 
quenched quark models we know of, thus making a 2 
charmonium resonance at 3.872 GeV very unlikely |14| . 
For further information and more references concerning 
A(3872), see e.g. a recent review [15], as well as our prior 
coupled-channel analysis |14j . 

The first suggestion of possible meson-meson molecules 
bound by pion exchange, in particular a DD* 1 state with 
quantum numbers 1 ++ or was due to Tornqvist |13| . 
With the discovery of A (3872) just below the D°D*° 
threshold, this idea was revived, of course. In the present 
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1 Henceforth, we omit the bar in DD , for notational simplicity. 
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paper, we intend to study the issue, not from Tornqvist's 
pion-exchange point of view, but rather as regards its 
possible implications for models based on quark degrees 
of freedom. In this context, it is worthwhile to quote from 
Ref. in which a molecular interpretation is advocated 
(also see Ref. |17|): 

"Independent of the original mechanism for 
the resonance, the strong coupling transforms 
the resonance into a bound state just below the 
two-particle threshold if a > or into a vir- 
tual state just above the two-particle thresh- 
old if a < 0. If a > 0, the bound state has a 
molecular structure, with the particles having 
a large mean separation of order a." 

(Note that, here, a is the 5*-wave scattering length.) Also: 

"In this case [1 ++ ], the measured mass Afj 
implies unambiguously that X must be either 
a charm meson molecule or a virtual state of 
charm mesons." 

In face of these peremptory claims about the molecular 
picture, we believe it is of utmost importance to study in 
detail the A (3872) wave function for a model in which the 
mechanism generating the meson is quark confinement 
combined with strong decay. Thus, we employ a sim- 
plified, coordinate-space version of the coupled-channel 
model used recently [14] to describe A (3872) as a uni- 
tarized and mass-shifted 2 3 P 1 charmonium state. The 
model's exact solvability then allows to obtain analytic 
expressions for the wave-function components, and fol- 
low bound-state as well as resonance poles in different 
Riemann sheets. 

The model is outlined in Se c. [TT] with details moved to 
Appendices [X] and [5] Section III is devoted to 5-matrix 
pole trajectories as a function of the two free parameters. 
In Sec. |IV| the two-component wave function is analyzed 
for several parameter sets, and is then used in Sec. [V|to 
compute cc probabilities and root-mean-square (r.m.s.) 
radii. The dynamical vs. confinement nature of the poles 
is discussed in Sec. lVH Conclusions are drawn in Sec. lVIII 



straightforward fashion. So for the purpose of the present 
paper, we resort to the equivalent |19| coordinate-space 
coupled-channel formalism of Ref. |5D] , which was used to 
study the influence of strong decay channels on hadronic 
spectra and wave functions, besides several more specific 
phenomenological applications. Furthermore, since here 
we do not aim to describe all aspects of the AT (3872) 
resonance in a detailed way, but rather want to focus 
on the importance of the charmonium component in the 
wave function, we restrict ourselves henceforth to a sim- 
ple two-channel system, viz. a 3 Pi cc state coupled to an 
S'-wave D°D*° channel, the latter being the dominant 

one, observed in the D°D ir° mode [5J. A partial study 
of this problem was already carried out in Ref. |21j . 

Consider now a system composed of a confined qq chan- 
nel coupled to a meson-meson channel M\M%. Confine- 
ment is described through a harmonic-oscilator (HO) po- 
tential with constant frequency to, having spectrum 



E= (2v + l c + -)w + m,+! 



(1) 



where v is the radial quantum number, l c the qq orbital 
angular momentum, and m q — rriq = 2fi q the constituent 
quark mass. In the scattering channel, no direct interac- 
tions between the two mesons are considered, with [if and 
If being the reduced two-meson mass and orbital angu- 
lar momenta in the free channel, respectively. Transitions 
between the two channels are modeled via an off-diagonal 
delta-shell potential with strength g, which mimics string 
breaking at a well-defined distance a. The correspond- 
ing Hamiltonian, transition potential, and 2x2 matrix 
Schrodinger equation are given in Eqs. pJ-((5]), with the 
usual definition u(r) — rR(r), where R(r) is the radial 
wave function. The exact solution to these equations is 
derived in Appendix |A"| 
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II. THE COUPLED cc-D D*° SYSTEM 

In our previous work on A (3872) [T3], we described 
the state as a unitarized radially excited 1 ++ charmo- 
nium resonance, by employing the resonance-spectrum- 
expansion (RSE) (THj formalism. The RSE description of 
A (3872) amounted to coupling the most relevant meson- 
meson channels allowed by the Okubo-Zweig-Izuka (OZI) 
rule to a spectrum of bare 1 ++ cc states, and also the 
OZI-suppressed but experimentally observed [B] p°J/4> 
and toJ/tjj channels. Thus, we found that unquenching 
shifts the bare 2 3 Pi state more than 100 MeV downwards, 
which settles as a very narrow resonance slightly below 
or on top of the D°D*° threshold. However, the RSE 
approach does not allow to obtain wave functions in a 



V = 



h c V 
V h f 



2[i c a 



S(r — a) 



= E 



»/ 



(4) 



(5) 



Once the 1x15 matrix (cf. Eq. ( A13 1) has been 



constructed from the wave function, possible bound or 
virtual states as well as resonances can be searched for. 
These occur at real or complex energies for which S blows 
up, or equivalently, when cot Si f (E) — i (cf. Eq. (A12)). 



Thus, unlike the pure HO spectrum, which is a real 
and discrete set of energies, the "unquenched" spectrum 
given by Eq. [5] includes complex energies, too, some of 



3 



which even have no obvious connection to the original, 
"bare" levels. The corresponding poles, called dynam- 
ical, are the result of attraction in the meson-meson 
channel generated by the interaction with intermediate 
quark-antiquark states. The light scalar mesons /o(500), 
if o (800), /o(980), and a (980) [6] are archetypes of such 
resonances |22| . On the other hand, we designate by 
confinement poles the ones that can be linked straight- 
forwardly to the bare states. Nevertheless, we shall show 
below — as we already demonstrated in previous work — 
that these two types of resonances are not so distinct af- 
ter all, since minor parameter changes may transform one 
kind into the other. But regardless of the type of pole, 
the corresponding radial quantum number in Eq. ^ is 
related to the energy of the coupled-channel system by 



v{E) 



E — 2m c 
2u 



l c + 3/2 



(6) 



Only in the uncoupled case, that is, for g = 0, one recov- 
ers the original HO spectrum, with v = 0, 1, 2, .... 

Now we apply the formalism to the coupled cc-D°D*° 
system. The cc channel is assumed to be in a 3 Pi state, 
i.e., with l c = 1, implying the D°D*° channel to have 
If = or 2. Nevertheless, we shall restrict ourselves here 
to the S'-wave channel only, which will be strongly domi- 
nant, especially near threshold. The fixed parameters are 
given in Table [I] where the meson masses are from the 
PDG |6J, while the HO frequency u and the constituent 
charm quark mass m c are as originally determined in 
Ref. [23] and left unaltered ever since. Thus, from Eq. 
we get the lowest two HO states at Eq = 3599 MeV and 
Ei = 3979 MeV, respectively. The former should give 
rise — after unquenching — to the 1 3 Pi charmonium 
state Xci(l-P) 0) with mass 3511 MeV, while the latter 
is the bare 2 3 Pi state, which cannot so easily be linked to 
resonances in the PDG tables, though both X(3940) and 
AT (3872) are possible candidates, in view of their mass 
and dominant DD* decay mode [B]. However, AT (3940) 
may just as well be the, so far unconfirmed, 2 x Pi (P ) 
state h c (2P) fi"4] . 

The two remaining parameters, viz. the string-breaking 
distance a and the global coupling g, have to be adjusted 
to the experimental data. Nevertheless, these parame- 
ters are not completely free, as they both have a clear 
physical interpretation, albeit of an empirical nature. 
Thus, a is the average interquark separation at which 
Pi quark-pair creation/annihilation is supposed to take 
place, while g is the overall coupling strength for such 
processes. Note that we do not assume a particular mi- 
croscopic model for string breaking inspired by QCD, like 
e.g. in a very recent paper [53] . Still, the values of a found 



TABLE I: Fixed parameters. 



in the present work are in rough agreement with our prior 
model findings, and even compatible |25_ with a lattice 
study of string breaking in QCD [26 . Concretely, we 
have been obtaining values of a in the range 1-4 GeV -1 
(0.2-0.8 fm), logically dependent on quark flavor, since 
the string-breaking distance will scale with the meson's 
size, being smallest for bottomonium. As for the cou- 
pling parameter g, its empirical value will depend on a. 
but also on the set of included decay channels. In realistic 
calculations, values of the order of 3 have been obtained 
(see e.g. our previous paper [T3] on the X(3872)). 



III. POLES 

The crucial test the present model must pass is its ca- 
pability of generating a pole near the D°D*° threshold. 
Indeed, a dynamical pole is found slightly below thresh- 
old for different combinations of the free parameters a 
and g, several of which are listed in Table [TTJ Examples 
are here given of bound states, virtual bound states, and 
below-threshold resonances, the latter ones only occur- 
ring for S'-wave thresholds as in our case. Note that poles 
of both virtual bound states and resonances lie on the 
second Ricmann sheet, i.e., the relative momentum has a 
negative inaginary part. From this table we also observe 
that larger and larger couplings are needed to generate 
a pole close to threshold when a approaches the value 
3.5 GeV -1 . We shall see below that this is due to the 
nodal structure of the bound-state wave function. 

Although a dynamical pole shows up near the D°D*° 
threshold, there still should be a confinement pole con- 
nected to the first radial 3 Pi excitation at 3979 MeV. 
Well, we do find such a pole, for each entry in Table [H] 
In Table |III| a few cases are collected, with the parame- 
ters tuned to generate a dynamical pole at precisely the 



TABLE II: Bound states (BS), virtual bound states (VBS), 
and resonances closest to threshold, for various g and a com- 
binations. 



Parameter 
Value (MeV) 



190 



15621 



m D o 
1864.8( 



2006.' 



m D o + m D ,o 
3871.84 



a (GeV" 1 ) 


9 


Pole (MeV) 


Type 


2.0 


1.149 


3871.84 


VBS 


2.5 


1.371 


3871.84 


VBS 


3.0 


2.142 


3871.84 


VBS 


3.1 


2.503 


3871.84 


VBS 


3.2 


2.531 


3871.84 - £12.01 


resonance 


3.3 


3.723 


3871.84 - i 4.45 


resonance 


3.4 


7.975 


3871.84 -i 0.39 


resonance 


3.5 


oo 






2.0 


1.152 


3871.84 


BS 


2.5 


1.373 


3871.84 


BS 


3.0 


2.145 


3871.84 


BS 


3.1 


2.507 


3871.84 


BS 


3.2 


3.083 


3871.84 


BS 


3.3 


4.194 


3871.84 


BS 


3.4 


8.254 


3871.84 


BS 


3.5 


oo 
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TABLE III: Pole doubling: pairs of poles (in MeV) for some 
sets of a and g values, chosen such that the dynamical pole 
settles at the X(3872) PDG mass. 



a (GeV" 1 ) 


9 


Dynamical Pole 


Confinement Pole 


2.0 


1.172 


3871.68 


4030.50 - il36.51 


2.5 


1.403 


3871.68 


4063.27 - il24.07 


3.0 


2.204 


3871.68 


4101.48 - i 88.03 


3.4 


8.623 


3871.68 


4185.85 - i 20.63 



X(3872) PDG jS] mass of 3871.68 MeV. Note, however, 
that the associated confinement pole is not necessarily 
of physical relevance, since at the corresponding energy 
several other strong decay channels are open, which no 
doubt will have a very considerable influence and possi- 
bly even change the nature of both poles. As a matter of 
fact, in our prior paper |14j . with all relevant two-meson 
channels included, the X(3872) resonance was found as a 
confinement pole, whereas dynamical poles were only en- 
countered very deep in the complex energy plane, without 
any observable effect at real energies. So here we show 
these results only to illustrate that pole doubling may 
occur when strongly coupling S'-wave thresholds are in- 
volved, as we have observed in the past in the case of e.g. 
the light scalar mesons [22 and the 15*0(2317) [27 . The 
issue of confinement vs. dynamical poles will be further 
studied in Sec. IVII 

In order to better understand the dynamics of the dif- 
ferent poles, we plot in Fig. [T] pole trajectories in the 
complex energy plane as a function of the coupling con- 
stant g, and for three diferent values of a. For vanishing 
g, the dynamical pole acquires a negative infinite imag- 
inary part and so disappears in the continuum, whereas 
the confinement pole moves to the real energy level of the 
bare 2 3 P 1 state, i.e., 3979 MeV. As g increases, and for 
both a = 2.0 GeV -1 and a = 3.0 GeV -1 , the dynamical 
pole moves to the real axis below threshold, becoming 
first a virtual bound state and then a genuine bound 
state. Note that, in the latter case, the real part twice 
attains the V(3872) mass even before the pole reaches 
the real axis, but the corresponding imaginary parts are 
much too large as compared with experiment [6] , so only 
the bound state can be considered physical. Finally, for 
a = 3.5 GeV -1 the pole does never reach the real axis, 
which would require an infinite coupling. For the other 
parameter sets listed in Table [TTJ we find intermediate 
situations. Another feature we can observe for all tra- 
jectories is an initial attraction and subsequent repulsion 
between the dynamical and the confinement poles. 



IV. WAVE FUNCTION 

Now we are in a position to study the V(3872) bound- 
state wave function in several situations. We choose 
two values for the string-breaking para met er, viz. a = 
2.0 GeV" 1 and a = 3.0 GeV" 1 . In Table ITvl five different 



FIG. 1: Pole trajectories of dynamical (left) and confinement 
(right) poles as a function of g, for a=2.0 GeV -1 (top), 3.0 
GeV -1 (middle), and 3.5 GeV -1 (bottom), respectively. In 
the last case, there is no bound state near threshold. Note: 
(i) poles in Table. |Hl| are here marked by *; (h) arrows along 
curves indicate increasing g. 
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binding energies (BEs) are chosen with respect to the 
D°D*° channel, including the PDG [5] value labeled by 
X . We have computed and normalized (see Appendix |A| 
the two-component radial wave function R(r) was com- 
puted and normalized for each of the five cases. In Fig. [2] 
we depict the cases labeled by A, X and D, the other two 
representing intermediate situations. A general feature 
we immediately observe is the typical S'-wave behavior 
of the D°D* a wave-function component i?/, while the cc 
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TABLE IV: Five chosen binding energies (BEs) in the D°D*° 
channel, for two different a values and the corresponding cou- 
plings g. 



TABLE V: Probabilities (in %) of the two wave-function 
componenents, for the cases specified in Table |fV| 



a 


(GeV- 1 ) 


2.0 


3.0 


Label 


BE (MeV) 


9 


Pole 


9 


Pole 


A 


0.00 


1.152 


3871.84 


2.145 


3871.84 


B 


0.10 


1.167 


3871.74 


2.191 


3871.74 


X 


0.16 


1.172 


3871.68 


2.204 


3871.68 


c 


1.00 


1.207 


3870.84 


2.311 


3870.83 


D 


10.00 


1.373 


3861.82 


2.899 


3861.84 



wave function R c is a P state, having one node. Also, 
\Rf\ is larger than |i? c | in most situations, for all r, except 
for unphysically large BEs (cf. plot D). Nevertheless, the 
two components are of comparable size for intermediate 
r values. Then, as the BE becomes smaller, the tail of 
Rf grows longer, as expected, whereas R c always be- 
comes negligible for distances larger than roughly 11-12 
GeV -1 . Now, the increased Rf tail affects the normaliza- 
tion of both R c and Rf. Thus, the ratio \Rf(r)\/\R c (r)\ is 
quite robust for most r values, as it does not significantly 
change with the BE. 



FIG. 2: Normalized two-component radial wave function R(r) 
for three BEs, corresponding to labels A, X, D in Table |IV| 
and two a values. Upper curves: Rf(r); lower curves: R c (r). 
Left: a = 2 GeV" 1 ; right: a = 3 GeV -1 . 
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V. PROBABILITIES AND R.M.S. RADII 

Having derived the AT (3872) wave function for several 
scenarios, we can now straightforwardly compute the rel- 
ative probabilities of the cc and D°D*° components (see 



a (GeV" 1 ) 


channel 


A 


B 


X 


C 


D 


2.0 
2.0 


cc 
D°D*° 


0.63 
99.37 


6.00 
94.00 


7.48 

92.52 


16.98 
83.02 


39.68 
60.32 


3.0 
3.0 


cc 
D°D*° 


0.97 
99.03 


9.01 
90.99 


11.18 

88.82 


24.65 
75.35 


55.54 
44.46 



Appendix [A| , with the results given in Table [vj for the 
five BEs and two a values from Table II VI Note that the 
probability in the D°D*° channel is only computed for 
normalization purposes, since in a more realistic calcu- 
lation at least the D ± D* T component would acquire a 
non-negligible probability as well, as the corresponding 
threshold lies only 8 MeV higher. Nevertheless, our sim- 
plification is unlikely to have an appreciable effect on the 
cc probability and will only increase that of the D°D*° 
component accordingly. Also note that the cc probabil- 
ity includes all 3 P± states, with the 2 3 Pi being domi- 
nant, because the corresponding bare eigenstate lies only 
100 MeV higher. However, also the 1 3 Pi state is non- 
negligible in the physical AT (3872) wave function. In the 
coupled-channel approach of Ref . [28J , a 1 3 P\ admixture 
of about 15% was found. Notice that — inevitably - 
unquenching not only mixes meson-meson components 
into the total bound-state wave function, but also quark- 
antiquark components of confinement states other than 
the one under consideration (also see Ref. Here, 
for a BE of 0.16 MeV, corresponding to the physical [6 
X(3872), case X in Table |Vj has a 7.48% cc probability 
for a = 2.0 GeV" 1 and 11.18% for a = 3.0 GeV" 1 . For 
smaller BEs, the cc probability decreases as expected, 
because of the growing weight of the D°D*° tail. On the 
other hand, for a BE of 10 MeV and a = 3.0 GeV" 1 , the 
charmonium probability becomes even larger than that of 
the meson-meson component. Now, the experimental er- 
rors in the average mass of the X(3872) and the D°D*° 
threshold allows for a maximum BE of 0.57 MeV, i.e., 
somewhere in between cases X and C. This would then 
correspond to a cc probability roughly midway in the 
range 7.48%-16.98% (a = 2.0 GeV" 1 ) or 11.18%-24.65% 
(a = 3.0 GeV -1 ). In the limiting case of zero binding, 
the cc probability would eventually vanish. Also notice 
that, in all five cases of Table |Vj the cc probability rises 
by about 50% when a is increased from 2.0 to 3.0 GeV -1 . 
Nevertheless, if we take a = 2.0 GeV -1 as in our Ref. |14j . 
we get a cc probability of 7.48%, very close the 7% found 
in Refs. [29l[3Q]. 



Next we use the normalized wave functions and 



Eq. (All) to compute the AT(3872) r.m.s. radius for the 
five cases discussed before (see Table IV I , with the results 
presented in Table VI It is interesting to observe that the 
r.m.s. radius, which in principle is an observable, is much 
less sensitive to the choice of a than de wave-function 
probabilities. Furthermore, the large to very large r.m.s. 
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TABLE VI: R.m.s. radii of the wave function, expressed in 
fm, for the cases specified in Table |TV| 



a (GeV -1 ) 


a (fm) | 


A 


B 


X C 


D 


2.0 


0.39 | 


100.22 


9.92 


7.82 3.10 


1.15 


3.0 


0.59 | 


100.14 


9.85 


7.76 3.05 


1.23 


radii in the 


various 


situations 


are 


hardly surprisin 





view of the small binding energies and the resulting very 
long tails of the D°D* a wave-function components (see 
Fig. [2] above). 

Using Eq. ( |A12l ), we now also evaluate the S'-wave scat- 
tering length 05 = — lim^^o [k(E) cot 60(E)] -1 . In case 
X and for a = 2.0 GeV -1 we thus find a, = 11.55 fm, 



which is large yet of the expected order of magnitude for 
a BE of 0.16 MeV. For even smaller BEs, the s catt ering 
length will further increase, roughly like oc l/\/BE. Let 
us here quote from Ref. |17| : 

"Low-energy universality implies that as the 
scattering length a increases, the probabilities 
for states other than D°D*° or D°D*° de- 
crease as 1/a . . . " 

Indeed, we verify from our Table [V] that — very roughly 
- the cc probability decreases as oc %/BE, and so like 
oal/a s . 



VI. STABILITY OF RESULTS AND NATURE 
OF POLES 

In this section we are going to study the stability of 
our results, as well as the nature of the found solutions. 
So let us vary the two usually fixed parameters, viz. u) 
and m c , in such a way that the bare 1 3 P± mass remains 
unaltered at 3599 MeV, whereas that of the 2 3 Pi changes 
as shown in Table. |VII| Thus, in case / Ei is lowered by 
25 MeV, while in case 77 it rises by 100 MeV. The trajec- 
tories for these two new situations are plotted in Fig. [3] 
For / we observe that, just as in the standard case de- 
picted in Fig. [T] two poles are found relatively close to 
the real axis, of a dynamical and a confinement origin, 
respectively. However, now it is the 2 3 Pi confinement 



TABLE VII: Probability of cc component and X(3872) r.m.s. 
radius for varying uj,m c , with bare Eq fixed at 3599 MeV, 
X(3872) pole at 3871.68 MeV, and a = 2.0 GeV" 1 . 





I 


standard 


II 


Ei (MeV) 


3954 


3979 


4079 


m c (MeV) 


1577.63 


1562 


1499.5 


cj (MeV) 


177.5 


190 


240 


g 


1.034 


1.172 


1.572 


cc (%) 


9.49 


7.48 


6.51 


Tr.m.s. (fm) 


7.72 


7.82 


8.83 



FIG. 3: Trajectories of dynamical and confinement poles. 
The bold curves represent cases I (top graph) and II (bot- 
tom graph) defined in Table |VII| and the others the stan- 
dard case of Fig. [TJ the solid (dashed) lines stand for normal 
(below-threshold) resonances. All trajectories lie on the sec- 
ond Riemann sheet. The pole positions for the g values in 
Table [VLT] are marked by *. 
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pole that moves steadily downwards and settles on the 
real axis below threshold, whereas the dynamical pole 
moves to higher energies and eventually approaches the 
real axis. So the poles interchange their roles when going 
from the standard case to case /. Nevertheless, the values 
of g needed to get a bound state at 3871.68 MeV are not 
very different in the two cases, viz. 1.172 vs. 1.034. Such 
a behavior was already observed almost a decade ago, 
namely for the Z?* (2317) [5] charmed-strange meson. In 
a first, two-channel model calculation [27 the £>* (2317) 
showed up as a dynamical resonance, settling below the 
S'-wave DK threshold, whereas the 1 3 Pq °s state turned 
out to move to higher energies, with a large width, simi- 
larly to the standard X(3872) case in Figs. [I] and [3] above. 
However, in a more complete, multichannel approach [31 
the situations got reversed, just as in the present case I. 
Also in our previous study |T3] of the X(3872), with nine 
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coupled channels, we reproduced the meson as a confine- 
ment pole. What appears to happen in the present case / 
is that shifting the bare 2 3 Pi state to somewhat lower en- 
ergies is just enough to make the confinement pole being 
deflected to the left and not the right when approaching 
the continuum pole. Clearly, there will be an intermedi- 
ate situation for which the left / right deflection will hinge 
upon only marginal changes in the parameters, but re- 
sulting in two completly different trajectories. Therefore, 
identifying one pole as dynamical and the other as linked 
to a confinement state is entirely arbitrary, the whole sys- 
tem being dynamical because of unquenching. At the end 
of the day, the only thing that really counts is where the 
poles end up for the final parameters. The trajectories 
themselves are not observable and only serve as an illus- 
tration how a coupled-channel model as the one empoyed 
here mimics the physical situation. Suffice it to say that 
the lower pole, representing the A(3872), is quite stable 
with respect to variations in the parameters, owing to its 
proximity to the only and most relevant OZI-allowed de- 
cay channel. The higher pole, on the other hand, should 
not be taken at face value, since a more realistic calcula- 
tion should include other important decay channels, such 
as D*D*, with threshold just above 4 GeV. 

Concerning the other scenario with changed parame- 
ters, labeled II in Table |VII| and depicted in the lower 
graph of Fig.[3j we see that the trajectories do not change 
qualitatively when going from the standard case to 77. 
There is a displacement of the right-hand branch, about 
100 MeV to the right on average, in accordance with the 
same shift of the bare 2 3 Pi state. But the change in the 
lower, dynamical branch, is much less significant, though 
the value of g needed to produce a bound state at 3871- 
68 MeV now increases to 1.572 (see Table |VII| . We also 
notice from Fig. [3] that the two pole-trajectory branches 
hardly move towards one another, signaling less attrac- 
tion between the poles due to a larger initial separation. 



FIG. 4: Normalized two-component radial wave function 
i?(r), for cases II and "standard", corresponding to param- 
eters in Table \VU\ Bold curves refer to case //, normal curve 
to R c for standard case. Note: Rf is indistinguishable within 
graphical accuracy for the two cases. 



Inspecting again Table VII as for the cc probability 
in cases I and II compared to the standard situation, 
we observe an increased value for case I and a decreased 
one for II. This is logical, since in case I the bare 2 3 Pi 
state lies closer to the X (3872), whereas in case II it lies 
farther away. Nevertheless, the difference in cc probabil- 
ity between / and II is only about 3%, i.e., less than the 
variation with a in the standard case. These comparisons 
lend further support to the stability of our results. 

Finally, in Fig. [4] we compare the wave function for 
case II with the standard one. We see there is no visible 
change in the Rf component. As for R c , the first maxi- 
mum gets somewhat reduced, but the secondary, negative 
bump even becomes a bit larger, owing to an inward shift 
of the node, lying now at about 3 GeV -1 . Yet, also in the 
case II the cc component is still very significant, despite 
the large separation of more than 200 MeV between the 
Y(3872) bound state and the bare 2 3 P i state. 

From the latter and all previous results we may safely 
conclude that the cc component of the X(3872) wave 
function remains nonncgligible in a variety of scenarios, 




being even of comparable size as the D°D*° component 
in the inner region, save at very short distances. 



VII. SUMMARY AND CONCLUSIONS 

In the present paper, we have employed a simple and 
solvable Schrodinger model to study the wave function of 
the A(3872) meson, by treating it as a coupled cc-D°D*° 
system with J PC = 1 ++ quantum numbers. Transi- 
tions between the two channels are described with the 
3 Po mechanism, through string breaking at a sharp dis- 
tance a. The exact solutions to the equations allow us to 
easily study the trajectories of 5-matrix poles as a func- 
tion of the decay coupling constant g. Thus, a dynami- 
cal pole is found, becoming a bound state just below the 
D°D*° threshold, for different string-breaking distances 
a, and an appropriate coupling g. On the other hand, 
the pole arising from the bare 2 3 Pi confinement state 
moves to higher energies and aquires a large imaginary 
part. However, the latter pole may not be very relevant 
physically, because of neglected additional meson-meson 
channels that will become important in that energy re- 
gion. 

As for the A (3872) radial wave function, the cc com- 
ponent R c turns out to be of significant size as compared 
to the D°D*° component Rf, especially for intermediate 
r values. Moreover, even for other trial BEs, the global 
shape of R c and its relative magnitude vis-a-vis Rf in 
the central region is remarkably stable. But the corre- 
sponding cc probability is relatively low, due to the very 
long tail of the D°D*° wave function at small binding. 
Quantitatively, for the average [6] Y(3872) binding of 
0.16 MeV, a cc probability of 7.5-11.2% is found, for a 
in the range 0.4-0.6 fm, which is compatible with other 
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recent approaches f29l 130] . The corresponding r.m.s. ra- 
dius turns out to be quite stable at about 7.8 fm, for 
the latter range of a values, while the iS-wave scattering 
length of 11.6 fm, for a ss 0.4 fm, is in agreement with 
expectations for a BE of 0.16 MeV. 

Finally, we have studied the nature of the 5-matrix 
pole giving rise to A (3872), by varying some of the oth- 
erwise fixed parameters. Thus, a drastic modification 
of pole trajectories is observed, for relatively small pa- 
rameter variations, making the A (3872) pole transform 
from a dynamical pole into one directly connected to the 
2 3 Pi bare confinement state. However, the correspond- 
ing changes in the cc probability and r.m.s. radius, as 
well as the coupling g needed to reproduce A(3872), are 
quite modest. 

In conclusion, we should revisit the claims about 
A (3872) made in Ref. [IS], quoted in the Introduction 
above, namely about the inevitability of A (3872) being 
a charm-meson molecule or virtual state, independently 
of the mechanism generating the state. Now, it is true 
that our analysis has confirmed some of the quantita- 
tive predictions in Ref. [16_, viz. concerning the vanish- 
ing probability of wave-function components other than 
D°D*° as the BE approaches zero, and the related be- 
havior of the D°D*° scattering length. However, we have 
also shown that the cc component is certainly not neg- 
ligible and quite stable, in a variety of scenarios. Es- 
pecially in electromagnetic processes, the prominence of 
this component for relatively small as well as intermedi- 
ate r values will no doubt result in a significant contribu- 
tion to the amplitudes. Moreover, as already mentioned 
above, the very unquenching of a 2 3 P\ cc state will not 
only introduce meson-meson components into the wave 
function, but also a contribution of the 1 3 Pi cc state, 
which can change predictions of electromagnetic transi- 
tion rates very considerably [25]. We intend to study 
such processes for A(3872) in future work, on the basis 
of a model as the one used in the present paper, by em- 
ploying the formalism developed and successfully applied 
in Ref. El]. 

For all these reasons, we do not consider A (3872) a 
charm-meson molecule, but rather a very strongly unita- 
rized charmonium state. As a matter of fact, we do not 
believe any non-exotic mesonic resonance — whatever 
its origin — qualifies as a true meson-meson molecule, 
simply because such a state will inexorably mix with the 
nearest qq states having the same quantum numbers. In- 
deed, we have demonstrated above that, even with a bare 
cc state 200 MeV higher in mass, the resulting cc compo- 
nent in the wave function is still appreciable. So let us 
conclude this discussion by quoting and fully endorsing 
the following statement from Ref. |16| : 

"Any model of the A(3872) that does not take 
into account its strong coupling to charm me- 
son scattering states should not be taken se- 
riously. " 



Acknowledgments 

One of us (G.R.) is indebted to R. M. Woloshyn 
for the invitation to a most stimulating miniworkshop 
at TRIUMF last year, where the topic leading to the 
present paper was debated. Thanks are also due to 
E. Braaten and K. K. Seth for very useful discussions 
on A (3872). This work was partially supported by the 
Fundagao para a Ciencia e a Tecnologia of the Ministerio 
da Educagao e Ciencia of Portugal, under contract no. 
CERN/FP/ 123576/20 1 1 . 



Appendix A: Solving the coupled-channel 
Schrodinger equation 

We twice integrate the Schrodinger equation ([5]) in or 
der to get two sets of boundary conditions, viz. Eqs. (IA1 



and (A2l 



u' c (r | a) - u' c (r I a) + ^"/( a ) = , 

u'Ar fa) - u'Ar | a) + ^u c (a) = 
1 ' aii c 



u c (r t a) = u c (r I a) , 
Uf(r f a) = Uf{r \. a) . 



(Al) 



(A2) 



A general solution to this problem is the two-component 
wave function given by Eqs. (A3 1 and (A4|, for the con- 



fined and meson- meson channel, respectively: 



u c {r) 



A c F c (r) r <a, 
B c G c (r) r> a; 



(A3) 



A f J lf (kr) , 
Uf(r) = { B f \j lf (kr)k 2l i +1 

cot (5i f (E)) - N lf (kr) 



r < a , 



r > a 



(A4) 

In Eq. (A3 1, the function F c (r) vanishes at the origin, 
whereas G c (r) falls off exponentially for r — > oo, their 
explicit expressions being 



F(r) = 



1 



r(Z + 3/2) 



J( (i+i)/2 e -«/2$(_ I/j i + 3/2, z ) , (A5) 



G(r) = ~T{- V )rz l l 2 e- Z l 2 ^!(~v, I + 3/2, z) , (A6) 

where $ and \1/ are the confluent hypergeometric func- 
tions of first and second kind (see Appendix [B| , respec- 
tively, rf— v) is the complex gamma function, v is given 
by Eq. Q, and z = [lujr 1 . Note that only in the case of 
integer v, i.e., for g = 0, do $ and ^> reduce to the usual 
Laguerre polynomials for the three-dimensional HO po- 



tential. Furthermore, the functions J and N in Eq. (A4| 
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are simple redefinitions of the standard spherical Bessel 
and Neumann functions, i.e., Ji(kr) = k~ l rji(kr) and 
N[(kr) = k l+1 rni{kr). From the boundary conditions 
(Al) and (A2|, as well as the wave-function expressions 
(A3) and (A4|, we get, with the definition k = ka, 



G' c (r)F c (a)~F c (a)G c (a) = ^ f («)F c (a)^ , 



J[ f (K)N lf («) - J lf {n)N[ f (k) = l£lj lt («)F c (a) 



a (X c 



(A7) 

)^. 
B t 



Using now the Wronskian relations 

W(F c (a),G c (a)) = F c (a)G' c {a) - F' c {a)G c {a) = 1 , 

(A8) 

W(N lf ( K ),J lf (n)) = N lf ( K )J[ f (K)-N[ f ( K )J lf (K) = -1 , 



and continuity of the wave function at r = a (cf. Eq. ([A2j) ) , 
we can solve for three of the four unknowns A c , B Cl Af, 
and Bf. Note that Eqs. (Al) and (A2) are not entirely 



linearly independent, so that solving all four constants is 
not possible. This is logical, as the overall wave-function 
normalization does not follow from the Schrodinger equa- 
tion. Expressing all in terms of A c then yields 



F c (a) 
G c (a) 



Af = 



Ji f {K)G c (a) 



Ar 



B f = g -^J lf {K)F c {a)A c 
a n c 



(A9) 



Note that, in order to obtain the D°D*° wave function 
in the outer region, we must substitute cot (5;, (-E)J = i 
in Eq. (A4| (also see below). Finally, the normalization 



constant M of the total wave function is determined by 
computing 



dr \u{r)\ - 
o Jo 



dr (u 2 c (r) + u)(r)) = N 2 . (A10) 



Then, we can also calculate the root-mean-square radius 
f — \J (r 2 ) of the two-component system by 



1 



/ drr 2 (u 2 c {r) + u 2 {r)) 



(All) 



As for the 5-matrix poles corresponding to resonances, 
bound states, or virtual bound states, cot (Si f (E)) can 
be solved from continuity of U/(r) at r — a in Eq. ( |A4| , 
resulting in the expression 



cot (S h {E))=- 



g^kjf f (K)F c (a)G c (a) 



1 -1 



with the 1x15 matrix simply given by 
cot (Si f {E)) + i 



cot (S h (E)) -i 



THj{k) 

(A12) 



(A13) 



Real or complex poles can then be searched for numer- 
ically, by using Newton's method to find the energies 
for which cot (6i f (E)^ — i, on the appropriate Riemann 
sheet. 

Appendix B: Special functions, numerical methods, 
and kinematics 



The confluent hypergeometric functions 4> and ^ intro- 
duced in Appendix |A| are defined in Ref. [33], Eqs. (6.1.1) 
and (6.5.7), respectively. Thus, the function $ is easily 
programmed as a rapidly converging power series, while 
the definition (6.5.7) of iff in terms of $ and the gamma 
function T then also allows straightforward computation, 
by employing Gauss's multiplication formula for r(— v) 
(see Ref. [M], Eq. (6.1.20)) so as to map the argument — v 
to lying well inside the unit circle in the complex plane, 
whereafter a very fast converging power-series expansion 
of l/r(-i/) (see Ref. pH], Eq. (6.1.34)) can be applied. 

The integrals for wave-function normalization and 
computation of r.m.s. radii are carried out by simple 
Gauss integration, choosing increasing numbers of points 
on a finite interval for the cc channel, and an infinite 
one for D°D*°. Note that, in the former case, the wave 
function falls off fast enough to allow convergence for a 
finite cutoff, whereas in the latter a suitable logarithmic 
mapping is used. In both cases though, because of the 
wave- function cusp at r = a and in order to avoid numer- 
ical instabilities, the domain of integration is split into 
two pieces, with up to 16 Gauss points in the inner re- 
gion and 64 in the outer one, thus resulting in a very high 
precision of the results. 

Although the A (3872) bound state can reasonably be 
considered a nonrelativistic system, we still use relativis- 
tic kinematics in the D°D*° channel, since parts of the 
resonance-pole trajectories involve relatively large (com- 
plex) momenta. For consistency, the same is done for 
all energies. The manifest unitarity of the S matrix is 
not affected by this choice. Thus, the relative D°D*° 
momentum reads 



HE) 



(Bl) 



where T and P are the threshold {m D *a +m D o) and pseu- 
dothreshold (m£)*o — mn°) energies, respectively. The 
corresponding relativistic reduced mass is defined as 



Idk 2 
^ E ^2dE 



E 



TPY 

eJ) 



(B2) 



Note that in the cc channel the reduced mass is defined 
in the usual way, i.e., fi c — m c /2, owing to the inher- 
ently nonrelativistic nature of the HO potential and the 
ensuing wave function. 
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